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Introduction: The ubiquity of accretion disks 

around pre-main sequence and young main sequence 
stars having the potential to form planetary systems is 
now well established [1], However, unknown is an ac- 
curate estimate of the fraction of single stars with 
disks that have produced planetary systems. Theoreti- 
cal models of particle aggregation show that if parti- 
cles can grow from submicron to mm to cm in size, 
then the formation of planetesimals is possible in the 
time before the disk dissipates [2], The problem re- 
mains to understand how grains condense from nebu- 
lar gases, and how relic interstellar grains survive and 
are modified by their transport in the disk. If grains 
are lofted above the disk photosphere by processes 
such as winds, turbulent convection, or changes in ver- 
tical structure, the evolution of dust can be investigat- 
ed by observing the properties of the small (< lpm) 
grains in the optically thin disk surface layer or atmo- 
sphere. 

In this work, we examine the thermal emission 
from three Herbig Ae/Be stars of similar spectral type: 
HD150193 (A IV), HD100546 (B9V), and HD179218 
(B9). The stellar ages of the objects are: 2X10 6 yrs, 
10X10 6 yrs, and 0.5 X10 6 yrs for HD150193, 
HD100546 and HD179218, respectively [3], We chose 
these objects for study because they are isolated (not 
belonging to any know star forming region), have little 
to no active accretion, possess possible processing cir- 
cumstellar disks, and have varying degrees of silicate 
crystallinity [4], We assemble and model SEDs for 
each system using a passive reprocessing circumstel- 
lar disk to constrain the amount of crystals and flaring 
of the disk in each of the systems. 

Observations: Spectrophotometry of HD150193, 
HD100546, and HD179218 was obtained using the 
NASA Ames Research Center High efficiency Faint 
Object Grating Spectrometer (HIFOGS) [5], The 
spectrum of HD179218 was obtained on the 2.34-m 
telescope at the Wyoming Infrared Observatory using 
32-bit FORTH telescope software. The spectra of 
HD 150 193 and HD 100546 were obtained on the 4-m 
Blanco telescope at the National Optical Astronomical 
Observatories Cerro Tololo Inter-American Observa- 
tory in Chile. 

We also obtained 1.2 — 23.0 pm IR broad- and nar- 
row-band photometric observations of HD179218 on 
2001 April 30.42 UT using the University of Minneso- 


ta Mount Lemmon Observing Facility 1.52 m tele- 
scope and the UM multi-filter bolometer [6], 

SED Assembly: Data from the European Space 
Agency Infrared Space Observatory (ISO) Short- 
Wavelength Spectrometer (SWS) and the Infrared As- 
tronomical Observatory(IRAS) are used to expand the 
wavelength coverage of the spectral energy distribu- 
tion (SED) of each of our three HAEBEs. All of the 
data sets for each object were scaled to the HIFOGS 
spectra. 

Disk Model: We use a mineralogically enhanced 
version of the radiative, hydrostatic models of passive- 
ly irradiated circumstellar disk [7] (hereafter CG97) 
that had been developed to include some discrete min- 
erals [8] (hereafter C01). The CG97 disk is composed 
of two parts: an optically thin surface layer and an op- 
tically thick interior. C01 updated the CG97 model by 
adding a range of grain sizes and incorporating a sim- 
ple mineralogy through the use of laboratory deter- 
mined optical constants, including water ice, amor- 
phous olivine, and metallic Fe. In this work, we use 
the updated model of C01 with our own enhanced fea- 
ture of adding crystalline olivine grains into the opti- 
cally thin surface layer. We also divide the disk into 
an inner and outer region in which we can indepen- 
dently vary the crystalline to amorphous silicate ratio 
to provide the best fit to the SED. 

Crystalline Silicates: It is difficult to model the 
thermal emission from crystalline silicate grains [9], In 
this work, we choose to calculate the crystalline sili- 
cate emission from ellipsoids by elongating the crys- 
tals along one of the three crystalline axes [10]. Based 
on the location of the resonance peaks in the HIFOGS 
and ISO SWS spectra of HD 150 193 and HD 1792 18, 
we use an axis ratio of 10:1:1 to compute the optical 
efficiencies of the crystals. The crystalline silicates 
are not coated with ice since to do so requires mixing 
theory which eliminates the distinctive crystalline res- 
onances. We compute the thermal emission from the 
crystals in the disk surface layer of each object. By di- 
viding the disk into two regions, (inner and outer) we 
can vary the crystalline to amorphous silicate ratio to 
fit the model to the SED. 

Results: Our model results are displayed in Fig. 1 . 
To produce the best-fit model to the SEDs, HD100546 
and HD179218 are both modeled with a disk 150 AU 
in radius, while HD 150 193 is modeled with a disk 5 
AU in radius. The best fit SEDs for HD 100546 uses a 
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photosphere to gas scale height (H/h) of 4, H/h = 3 for 
HD150193, and H/h = 1 for HD179218. This means 
at a disk radius of 5 AU the disk around HD 100546 
flares about 80% more than the disk around 
HD179218, and about 25% more than the disk around 
HD 150 193. The model fits suggests that the measured 
emission from HD 150 193 is dominated by warm dust 
in a flaring disk close to the central star, and that the 
emission from grains in an extended disk (> 5 AU) is 
relatively small compared to the other two objects. 

HD100546 and HD179218 show evidence of crys- 
talline silicates, while the emission from HD 150 193 is 
dominated by amorphous silicates. For the two ob- 
jects with crystalline silicates, we find that a better fit 
is produced to the observed SEDs of the objects if we 
use a model in which a higher ratio of crystalline to 
amorphous silicates is located in the inner regions of 
the disk (< 5 AU) compared to the outer regions of the 
disk (5 - 150 AU). The inner region of HD100546 
has 30% more crystalline silicates compared to the 
outer region. This is contrary to the findings of [11] 
who used a spherical shell model to calculate a factor 
of almost 10 higher fraction of crystalline silicates in 
regions greater than 10 AU. The inner region of HD 
179218 has 76% more crystalline silicates compared 
to the outer region. 

All three SEDs are best fit using a grain size distri- 
bution with a slope of q = — 3.5 (i.e., a -35 ) . 

Discussion: It is difficult to make any statistically 
significant conclusions about disk evolution from 
modeling only three objects. However, we can make 
some interesting observations based on our results. 
Although there is no evidence for grain growth (all 
three SEDs are modeled with the same grain size dis- 
tribution) there is a difference in the radial distribution 
of crystalline silicates. A possible evolutionary sce- 
nario is one in which when circumstellar disks form 
around young stars, they are thought to be primarly 
composed of amorphous ISM grains [12]. As the disk 
evolves, crystalline silicates are condensed [13] or an- 
nealed [14] in the inner radial regions of the disk, ei- 
ther through heating [15] or through shocks in the disk 
[16]. The crystals are then transported to the outer re- 
gions of the disk [17]. Therefore, from this scenario, 
we can conclude that the disk around HD 150 193 is the 
least evolved since it does not contain any crystals, 
followed by HD 1792 18 which has an even larger 
amount of crystals in the inner radial region compared 
to the outer radial region, and finally HD 100546 
which has a larger amount of crystals in the outer re- 
gions of its disk compared to HD 1792 18. 
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Figure 1. C01 model computed for HD150193, 
HD100546 and HD179218. The disk interior (dotted 
line), disk surface (dashed line), and stellar blackbody 
(dot-dash line) are co-added to produce the model 
SED (solid line). The model SED is compared to the 
assembled data sets including: ISO SWS spectra (sol- 
id, jagged line), HIFOGS spectra (open circles) and 
IRAS photometry points (gray squares). HD 179218 
also has MLOF photometry points (filled circles). 




